Theory of Thermal Conductivity in Y Ba 2 Cu 3 07s 
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We calculate the electronic thermal conductivity in a d-wave superconductor, including both the 
effect of impurity scattering and inelastic scattering by antiferromagnetic spin fluctuations. We 
analyze existing experiments, particularly with regard to the question of the relative importance of 
electronic and phononic contributions to the heat current, and to the influence of disorder on low- 
temperature properties. We find that phonons dominate heat transport near T c , but that electrons 
are responsible for most of the peak observed in clean samples, in agreement with a recent analysis 
of Krishana et al. In agreement with recent data on YBa2(Cui- x Zn x )s07-s , the peak position is 
found to vary nonmonotonically with disorder. 



PACS Numbers: 74.72.Bk,74.25.Fy,74.62.Dh 

Introduction. Analysis of transport experiments in the 
superconducting state of the high- temperature cuprate 
superconductors has already provided the most com- 
pelling evidence for electronic pairing in these materials. 
The collapse of the quasiparticle relaxation rate below T c , 
as observed in optical |l],^j and microwave || measure- 
ments, is not observed in classic superconductors, and is 
most naturally interpreted in terms of a gapping of the 
spectral density of electronic excitations responsible for 
inelastic scattering just above T c . This collapse is now 
understood to be responsible, e.g., for the peak at inter- 
mediate temperatures in the microwave conductivity of 
YBa 2 Cu 3 O e . 95 . | 

Thermal conductivity measurements provide informa- 
tion on order parameter symmetry and quasiparticle re- 
laxation, and have the advantage that they are bulk 
probes not subject to extrinsic surface effects which have 
hampered the interpretation of the low-T microwave con- 
ductivity. H They have the disadvantage that the elec- 
tronic contribution to the heat current must be sepa- 
rated from the phononic one. As pointed out by Yu et 
al., Q the similarity between the microwave conductiv- 
ity peak and measurements of the thermal conductiv- 
ity k(T) in YBa2Cu30js single crystals suggests that 
at least part of the thermal conductivity peak should 
be due to the electronic thermal conductivity K e i(T), in 
contrast to earlier analyses of this peak in terms of a 
phonon conductivity K p h (T) alone. |g§ In the work of 
Yu et al., the phononic mean free path £ p h is assumed 
to vary only weakly with T < T c , whereas in the Pea- 
cor et al. approach, t v h is assumed to be dominated by 
phonon- electron relaxation, leading via the pair correla- 
tions in the electronic system to an exponential behavior 
below T c . In Ref. j5j it is assumed that Kph S> k e i over 
the entire temperature range, whereas Yu et al. deduce 
K P h{T c ) — 2 — 3n e i(T c ) using the measured <7i(T c ) on sim- 
ilar quality crystals, and assuming the Wiedemann- Franz 
law K e i = LqTo\ with the free electron Lorenz number 

u. 

In this paper, we adopt a theoretical model of elec- 
tronic transport in a d-wave superconductor limited 
by impurity and spin fluctuation scattering which has 



proven successful in describing many of the systematics 
of microwave measurements, and apply it to calculate 
the electronic thermal conductivity. We analyze exper- 
iments on Zn- doped Y ' Ba-iGu^O^ -$ to argue that a) 
phonons do in fact dominate heat conduction at T c ; b) 
a peak in n p h does indeed occur at about 20-25K; c) 
electronic conduction does nevertheless provide most of 
the peak in clean samples; and d) the temperature at 
which the peak in k(T) occurs may vary nonmonotoni- 
cally with disorder. Conclusions a)-c) agree with a recent 
thermal Hall conduction measurement and analysis by 
Krishana et al. jjj Conclusion d) agrees well with recent 
data on Y Ba2{Cui- x Zn x ) 3 0-!-s over a wide range of 
Zn concentrations. || We further focus on the very low- 
temperature behavior of the thermal conductivity in the 
d-wave model, discussed recently in considerable detail 
by Graf et al. || In particular, we analyze the disorder 
and phase shift dependence of the "universal" linear-T 
term in n predicted recently by several groups. JlO|,p 11 
Since at very low temperatures the phonon mean free 
path has saturated, this contribution is a direct reflec- 
tion of electronic correlations and relaxation. 

Electronic thermal conductivity. The electronic ther- 
mal conductivity k for an unconventional superconduc- 
tor (lj [H| is evaluated using a Kubo formula for the 
heat-current response as in the original treatment for an 
s-wave superconductor. [fl6f In this letter we particularly 
wish to study the d x 2_ y 2 pair state thought to provide 
a good description of optimally doped Y BaiCu^Oi-b 
JlT| , p~8[ | . For simplicity we work in what follows with the 
approximate order parameter — Ao cos(2</>) over a cir- 
cular Fermi surface to describe ab plane transport. The 
impurity-averaged matrix (Nambu) electron propagator 
in such a state is given by 



5(k,w) 



CkZ 3 + iA^t 1 
-€-|A k | 2 



(1) 



where a_ l and r l are the Pauli matrices in spin and 
particle-hole space, respectively. Here, we have already 
exploited the assumed particle-hole symmetry of the nor- 
mal state, as well as the symmetries of the gap functions 
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which lead to vanishing renormalizations for both the or- 
der parameter and the single-particle energies. In this 
case, only self-energy contibutions to the frequency u>, 
namely Q = to — So need to be included. p3 The self- 
energy Eo due to the elastic impurity scattering is treated 
in a self-consistent i-matrix approximation and is given 
by E = TG Q /(c 2 - Go), where T = mn/{-KN ) is the 
normal state unitarity limit scattering rate depending on 
the concentration of defects n^, the electron density n, 
and the density of states at the Fermi level Nq. The 
quantity c = cot 8$ parameterizes the scattering strength 
of an individual impurity through the s-wave phase shift 
6q. In this work we consider only near-unitarity limit 
scattering c ~ since it is clear that weak scattering 
will lead to a weak temperature dependence inconsistent 
with experiment for the states in question. The inte- 
grated propagator is Go = (l/2TrN )^2 k Tr{T°g(k,ui)}. 
The equation for the self-energies are then solved self- 
consistently together with the gap equation A^ = 
-TE U „ £k< V kk -(l/2)Tr{r 1 g(k>„)} . 

The above approximation is insufficient to describe 
transport at temperatures close to T c , where inelastic 
scattering is known to dominate. As in Refs. p9| , p0| , we 
adopt a model of scattering by antiferromagnetic spin 
fluctuations based on an RPA treatment of the Hubbard 
model with parameters chosen to reproduce normal state 
NMR and resistivity data in YBa 2 Cu 3 7 s. |HJ The re- 
laxation rate due to spin fluctuations 1 / Ti n is quasilincar 
in temperature above T c , and falls as ~ T 3 in the su- 
perconducting state due to a) a crossover to Fermi liq- 
uid behavior below a spin fluctuation scale, and b), an 
additional factor of T due to the restriction of relevant 
quasiparticle momenta to the vicinity of the d-wave or- 
der parameter nodes. (2l]] To include inelastic scattering 
in the model in a crude way, we make the replacement 
Eo —> E — i/2r in . Although we have adopted a particular 
microscopic model, we emphasize that the two features 
a) and b) noted may be common to many models of in- 
elastic relaxation, which will then lead to qualitatively 
similar results. 

The bare heat current response is now given by a con- 
volution of the Green's function g with itself at zero ex- 
ternal frequency and wave vector weighted with the bare 
heat current vertex ivvp klf.- fl6| Impurity scattering ver- 
tex corrections to current-current correlation functions at 
q = have been shown to vanish identically for even par- 
ity states (Ak = A_k). |l4| For the diagonal thermal 
conductivity tensor one obtains 



<i(T)/T 



G 



dui ( — 



dui 



Ki(u,T) (2) 



^,T) = ^Re(k--" 2 + l"l 2 - 2 l Ak ' 2 



, (3) 



where &' and u>" are the real and imaginary parts of u>, 
f is the Fermi function, and T to t = T + l/2ri„ is the 



total quasiparticle scattering rate. We have numerically 
evaluated Eqs. (2-3), and show results in Figure 1. In 
the clean limit, the combination of the relaxation rate 
collapsing with decreasing temperature due to gapping 
of the spin fluctuation spectral density and the rapidly 
decreasing number of quasiparticles at low T leads to 
a peak in the thermal conductivity very similar to that 
found for the electrical conductivity. As impurities 

are added, the collapse of the inelastic scattering rate is 
cut off at progressively higher and higher energy scales, 
such that the peak moves to higher temperatures and 
simultaneously weakens. A preliminary report on these 
findings is contained in Ref. p4| . 
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FIG. 1. Normalized electronic thermal conductivity 
k/kjv for normalized impurity scattering rate T/T cQ 
= 0.001,0.005,0.01,0.05, and 0.25, A /T c = 3 and 
l/T in (T c ) = T c 

Phonon thermal conductivity. Were the phonon and 
electron thermal conductivities comparable at T c , the 
electronic peak in clean samples would, according to Fig- 
ure 1, be very large, leading to ratios K{T pea k) I k(T c ) 
much larger than the observed range of about 1.8-2.4. 
P,^|p4|j8]| We therefore expect on theoretical grounds that 
the phononic conductivity at T c is several times larger 
than the electronic conductivity. (2^] This conclusion was 
reached independently by Krishana et al., [Q who per- 
formed a semiclassical Boltzmann type analysis of their 
measurements of the diagonal and off-diagonal compo- 
nents of the tensor K%j in a magnetic field, assuming as 
the origin of the off-diagonal terms skew-scattering of 
quasiparticles off the vortex lattice. The diagonal phonon 
conductivity K p h was then extracted from the data under 
the assumption that only electrons were skew scattered. 
Here we adopt the n p h as shown in Figure 2 following the 
analysis of Ref. , using a value of the transverse scatter- 
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ing cross section such that n{T pea k) — k(T c ) is primarily 
of electronic origin. JjJ We use the form of n p h shown for 
all calculations at any impurity concentrations, thereby 
neglecting the weak effect of point defects on the long 
wavelength phonon mean free path, T~ oint ~ w 4 . The 
overall scale of n p h is set by its value at T c , which is 
found to be roughly seven times n e i(T c ). 
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FIG. 2. Solid line: electronic thermal conductivity re e ;; 
dash-dotted line: phononic thermal conductivity k v h\ dashed 
line: K vs. T/T c for T/T c0 = 0.007. Data are from Ref. 7. 

Zn substitution: high temperatures. In Figure 2, we 
plot the theoretical k(T) obtained by adding k p h(T), de- 
termined as above, to K e i(T), with disorder parameter 
r chosen to give rough agreement with the peak height 
and position of the nominally clean sample of Ref. M. 
We note that the quasiparticle mean free path extracted 
from nominally pure crystal data in Ref. [Q, is actually 
similar to the mean free path in the 0.15% Zn sample of 
Ref. ||, making the assignment T/T c — 0.007 consistent 
with earlier analysis. P0[] In Figure 3 we now show re- 
sults for the total thermal conductivity with systematic 
Zn doping compared to data of Ref. (^] on twinned crys- 
tals of Y Ba,2(Cu\- x Zn x )iO-!-b. Note that the impurity 
scattering parameters T for the various curves were cho- 
sen to reflect Zn concentrations of 0.1%,0.2%,0.7%,1.7% 
using the identification of a contribution to the scattering 
rate T/T c q = 0.12 per 1% Zn. This is a factor of 2 larger 
than the contribution extracted from comparisons with 
microwave data in Ref. pcfl , but consistent with the error 
bars cited in that work. We have also added in all cases 
a residual scattering rate of T/T c = 0.007 consistent with 
the nominally pure sample as in Figure 2, apparently rep- 
resentative of oxygen defects in the near-optimally doped 
samples. With this assignment, there are no further free 
parameters in the theory. We note first that the posi- 
tion of the peak in temperature is nonmonotonic in both 
the data and the theory. In the current theory this occurs 



because, while the electronic peak initially dominates and 
moves upward in T with disorder, the phonon peak at 20- 
25K eventually becomes more important as the scale of 
K e i is reduced by disorder. We expect that continued Zn 
doping will lead to saturation and a peak position fixed 
at 20 — 25-fT when n e i disappears, leaving the (roughly 
impurity independent) K p h of Fig. 2. 
O 
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FIG. 3. Total thermal conductivity k vs. T[K]. Sym- 
bols are data from Ref. 8, Zn concentration 0.0%, 
0.1%,0.2%,0.7%,1.7%. Solid lines: theoretical k for 
r/Tco = 0.007 + 0.12(%Zn). 

Wiedemann- Franz and anisotropy ratios. Measured 
values of the ratio k(T c ) /T c a(T c ) are typically 1.0—1.5 x 
10~ 7 IT - VL/K 2 . If we take R = K ph (T c ) / n el (T c ) = 
7 as above, we find for the Lorenz number L = 
K e i(T c )/Ta(T c ) = 1.2 - 2.0 x 10~ 8 IT - fl/K 2 . Given 
that electron correlations can change this ratio substan- 
tially, this is quite close to the free electron value of 
L = 2.44 x 10~ 8 W - fl/K 2 . The measured anisotropy 
ratio in the plane, Kb(T c ) / K a (T c ) is 1.2-1.3 for clean un- 
twinned crystals. If we assume that all anisotropy 
arises from the electronic component, and using R ~ 7 
once again, we find Hb,ei(T c ) / K a ,ei(T c ) = 2.6 — 3.4. This 
is close to the measured electrical conductivity ratio of 
Ob jo a — 2.4 reported by Zhang et al. j2^|. 

Zn substitution; low temperatures. Direct information 
on electronic properties may be obtained by working at 
very low temperatures, such that the phonon contribu- 
tion, which should fall as T 3 when the mean free path 
saturates due to boundary scattering, is negligible. In 
the resonant scattering limit considered here, a signif- 
icant linear-T electronic contribution n e i ~ aT should 
dominate the phonon conductivity in a d-wave supercon- 
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ductor. The thermal conductivity in this limit is "uni- 
versal" in the sense that the prefactor a is independent 
of the impurity scattering rate to leading order, ]To|.p|,pl| 
in analogy to the limiting T = 0,uj —> conductivity 
coo = (ne 2 )/(rmrA ). p3] As pointed out by Graf et al. 
the Wiedemann-Franz law is obeyed exactly at T = in 
the clean limit, i.e. a — > Lqo-qq. 
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FIG. 4. Normalized electronic thermal conductivity 
K e i/ LoaooTvs.T /T c for varying electronic scattering phase 
shifts, c = 0,0.1,0.13, r/T c = 0.007. Insert: K el /L a 00 T 
vs. r/T c for fixed T = 0.01T C for c = (solid line); c = 0.1 
(dashed line); and c = 0.3 (dashed-dotted line). 

While coo is difficult to measure due to its small size 
and to possible extrinsic contributions alluded to above, 
the linear term in k(T) should be clearly visible, and its 
magnitude and impurity dependence sensitive tests of the 
"dirty d-wave" model. 

The exact expression for the limiting value of K e i/T 
fo r T < T c is fny§ a = L <T 0Q k''E(k'), where k' = 
\J Ao/(Aq + 7 2 ), E is a complete elliptic integral of the 
2nd kind, and 7 = — ImSo(w = 0). In the unitar- 
ity limit, c = 0, the K e i = aT relation holds over a 

temperature range T^^/TAq, whereas away from this 
point the range of validity quickly vanishes. A quasi- 
linear behavior may nonetheless be observed; for ex- 
ample, in the opposite limit limit, c > 1, we have 
Ka(T)/T ~ [K e/ (T c )/T c ](2rr(T c ))- 1 above an exponen- 
tially small crossover scale ~ A exp — A (l + c 2 )/T. To 
illustrate the possible range of behavior we plot in Fig. 
4 Kei/T at T/T c = 0.01 for various values of c close to 
the unitarity limit, and as a function of Zn concentra- 
tion. The main point we wish to make here is that it 
is possible that small deviations from the unitarity limit 
may lead at low T to behavior quite different from that 
predicted for unitarity limit scattering. 



Conclusions. We have shown that "dirty d-wave the- 
ory" provides a good account of the systematic behavior 
of the thermal conductivity in Y Ba2(Cu\- x Znx)30-r-s- 
In particular, it allows for a simple understanding of the 
size and position of the intermediate temperature peak 
as a function of disorder. The analysis is consistent with 
earlier studies of the microwave conductivity, enabling 
semiquantitative predictions. On this basis we have pro- 
vided strong evidence for a phononic conductivity signifi- 
cantly larger than its electronic counterpart near T c . The 
electronic component is nevertheless responsible for the 
peak in clean samples, but disappears with per cent level 
Zn doping. Finally, we have discussed the low-T limit- 
ing behavior of n, and pointed out that small deviations 
from the unitarity scattering limit can give rise to quite 
nonuniversal results. 
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